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We report the first observation of the decay B"*" ip{3770)K'^ where '0(3770) is reconstructed in 
the D°D° and D+D' decay channels. The obtained branching fraction is B{B+ V(3770)/S'+) = 
(0.48±0.11±0.07) X 10"^. We have measured the branching fraction for the decay B+ D°D°K+ 
to be (1.17 ± 0.21 ± 0.15) x 10"^. We have set a 90% confidence level upper limit of 0.90 x 10"^ 
for the decay B'^ — » D~ . We also present the results of a search for possible decays to 
DD and D^ZJ^tt" of the recently observed X(3872) particle. The analysis is based on 88fb~^ of 
data collected at the T(4S') resonance by the Belle detector at the KEKB asymmetric-energy e"'"e~ 
collider. 



PACS numbers: 13.25.Hw, 13.20.Gd 

B decay modes with charmonium in the final state 
are extensively used by the Belle and BaBar Collabo- 
rations for measurements of the CP violation parameter 
sin2(^i 0,0 . Belle has recently reported the first observa- 
tions of the decays XcoK'^ j3(| and B Xc2X 0. 
The decay rates for these modes were measured to be 
comparable to those for J/^ and tp{2S). 

In contrast to the charmonia seen so far in B decays, 
the ^(3770) state is just above open charm threshold 
and decays dominantly to pairs of D mesons |^. The 
'0(3770) is generally considered to be predominantly the 
l^Di charmonium state. However, it has a non-zero lep- 
tonic width, which indicates that there is some mixing 
with the nearby ip{2S) S-wave state 0. A large S-D-wave 
mixing angle could result in comparable decay rates for B 
decays to the i/'(3770) and the ^{2S). For a pure D-wave 
state, an estimate of B{B '!/)(3770)X) based on the 
color-octet model gives a value of 0.28% 0, which is as 
large as the measured values for J/'0 and tp{2S). Exper- 
imental studies of '0(3770) production in B decays test 
theoretical models and provide additional information on 
the structure of the 0(3770) wave-function. 

In this paper, we report the first observation of the de- 
cay B"*" — > 0(3770) if"*" 0. We also report measurements 
of the B+ D"D°K+ and B+ D+ D' K+ decay 
modes lil and searches for the decays B+ X(3872)fs:+, 
X(3872) -> DD {D°D^iT°). The analysis is performed 
using data collected with the Belle detector at the KEKB 
asymmetric-energy e^e~ collider |lO|. The data sample 
consists of 88 fb""'^ taken at the T(45) resonance, which 
corresponds to 96 x 10^ BB pairs. 

The Belle detector is a general-purpose spectrometer 



with a 1.5-T superconducting solenoid. Charged par- 
ticle tracking is performed by a silicon vertex detector 
(SVD) composed of three concentric layers of double 
sided silicon strip detectors, and a 50-layer drift cham- 
ber (CDC). Particle identification for charged hadrons 
is based on the combination of energy-loss measurements 
(dE/dx) in the CDC, time-of-flight measurements (TOF) 
and aerogel Cherenkov counter (ACC) information. A 
CsI(Tl) electromagnetic calorimeter (ECL) located inside 
the solenoid coil is used for the detection and identifica- 
tion of photons and electrons. The detection system for 
K1 mesons and muons consists of alternating layers of re- 
sistive plate counters and 4.7-cm-thick steel plates. The 
Belle detector is described in detail elsewhere [llj. 

We select charged pions and kaons that originate from 
the region dr < 1 cm, \dz\ < 3 cm, where dr and dz are 
the distances of closest approach to the interaction point 
in the plane perpendicular to the beam axis and along the 
beam direction, respectively. Charged kaons are required 
to satisfy C{K) / {C{K) + C{Tr)) > 0.6, where C{K/Tr) is 
the particle identification likelihood for the K/tt hypothe- 
ses calculated by combining information from the TOF, 
ACC and dE/dx measurements in the CDC. Candidate 
tt'^ mesons are identified as pairs of non-charged-track- 
associated ECL clusters that have an invariant mass 
within ±15 MeV/c^ of the tt" mass. To reduce combina- 
torial background, the energy of each photon is required 
to be greater than 50 MeV and the momentum of the 
TT*^ in the center of mass system (cms) is required to be 
greater than 0.15 GeV/c. 

The D^ meson is reconstructed in the K^n^, 
if"7r+7r+7r^ and X^tt+tt^ modes, and the D^ in the the 
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TABLE L Summary of the fit results, ofRciencies, statistical significance and branching fractions for D D and 

B+ D+D~K+ decays. 



Mode 


AE yield 


Efiiciency(10~'') 


S(10"^) 


Significance 




97.5 ± 17.6 


8.7 


1.17 ±0.21 ±0.15 


5.5 a 


B+ ->■ D+D-K+ 


20.7 ±9.9 


5.0 


< 0.90 (90% CL)(0.43 ± 0.21 ± 0.06) 


2.7 a 



K 7r+7r+ and K 7r+ modes. We use a ±10MeV/c^ 
D signal window for the charged modes 2.5 cr) and 
±15MeV/c2 for K-tt+ti° mode (~ 2 a). Mass- and 
vertex-constrained fits are apphed to aU D candidates 
to improve their momentum resolution. The candi- 
dates {i.e. DD pairs combined with the positive kaons in 
the event) are identified by their cms energy difference, 
IS.E — TiiEi — -Ebeamj and their beam-energy constrained 
mass, Afbc = v/^^b^^Ii^i^^T^i^iF: where E'boam = %/s/2 
is the beam energy in the cms and pi and Ei are the 
three-momenta and energies of the i?+ candidate's de- 
cay products. For the D'^D'^K^ final state, we 
require that one is reconstructed in the Z?" K~tt~^ 
mode, which has the smallest background. We accept B 
candidates with 5.272 GeV/c^ < M^c < 5.288 GeV/c^ 
and |Ai5| < 0.2 GeV. To suppress the continuum back- 
ground we require the ratio of the second to the zeroth 
Fox- Wolfram moments R2 to be less than 0.5 and 
|cos6'thr.| < 0.8, where ^thr. is the angle between the 
thrust axis of the B candidate and the thrust axis of 
the rest of the event. The last cut is not applied for the 
cleanest subset of B candidates where both Z?°'s are re- 
constructed in the Ktt mode. In the case of multiple B 
candidates, we choose the candidate with the smallest 
value of = ((^bc - MB+)/(TMbJ^- 

The AE distributions for the B+ D^D^K+ and 
B+ D^D^K^ candidates are shown in Fig.Q] where 
the superimposed curves are the results of the fits. The 
fit to the AE distribution is a sum of a Gaussian with a 
fixed width taken from Monte Carlo (MC) to describe the 
signal and a first order polynomial to parameterize the 
background the fit to the AE distribution, the 

region AE < —0.08 GeV is excluded to avoid contribu- 
tions from other B D^*^D^*^K decays. Tablenjsum- 
marizes the results of the fits, the reconstruction efficien- 
cies , statistical significance of the signals and the 
calculated branching fractions. For the latter, we assume 
N{B°B") = N{B+B-). For the D+D-K+ final state, a 
substantial signal is not seen and we set a 90% confidence 
upper limit. The systematic error in the branching frac- 
tion measurement is dominated by the uncertainty in the 
tracking efficiency (1% per track), kaon identification ef- 
ficiency (2% for each kaon), tt" reconstruction efficiency 
(6%), branching fraction uncertainty (in total 8%), 
MC statistics (3%), the signal and background parame- 
terization (5%). 

We plot the D^D^ and D^D^ invariant mass dis- 
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FIG. 1: The AE distributions for the B+ 
D° D° K+ {u-pper) and B+ D+ D' K+ {lower) candidates. 
Points with errors represent the data and curves show the 
results of the fits described in the text. 



tributions for events in the B signal region defined as 
5.272 GeV/c^ < Khc < 5.288 GeV/c^ and |A£;| < 
0.02 GeV in F_ig. I^a) and Fig.l^Jb), respectively Here, 
for D^D^K^ candidates, when one of the IJ^'s is 

reconstructed in the K^h^tt^ mode, we use a looser AE 
cut (|Ai?| < 0.025 GeV) to take into account the poorer 
energy resolution due to shower leakage in the ECL. In 
the case of multiple B candidates, we choose the can- 
didate with the smallest value of = {AE/a^^E)"^ ± 
((Afbc - Mb+)I(Tm^,)^- The M{D°D°) distribution has 
a peak at low masses, which we attribute to the 7/; (3770) 
signal. 

The superimposed hatched histogram in Fig. [2ta) 
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TABLE II: Result of the measurement of Am = M(i/'(3770)) — Ad{tp{2S)) obtained in tiiis paper and previous measurements. 



Belle MARK I DELCO MARK II 

Am = M(V'(3770)) - M{iP{2S)), MeV/c^ 92.4 ± 3.0 ± 1.3 88 ± 3 86 ± 2 80 ± 2 



shows the M{D'^D^) mass distribution for events in the 
/S.E sidebands The curve (shown as solid line) is 

the result of a fit where the low-mass peak is described 
by a relativistic p-wave Breit-Wigner function ^| with 
a floating mass and its natural width fixed to its nom- 
inal value of r(V^(3770)) = 23.6 MeV/c^ Q. The com- 
binatorial background along with the contribution from 
non-resonant 5+ D^D^K'^ decays is described by the 
product of left and right two-body threshold functions 
and is represented as a dashed line in Fig. |2Ia). 

The fit yields a -0(3770) signal of iV = 33.6 ± 8.3 
events with a statistical significance of 5.9(t. The mass 
of the 1/1(3770) is found to be Af (0(3770)) = 3778.4 ± 
3.0 ± 1.3 MeV/c^, which corresponds to a mass differ- 
ence Am = M(0(377O)) - M(0(25)) = 92.4 ± 3.0 ± 
1.3 MeV/c2, where we used M(0(25)) = 3685.96 ± 
0.09 MeV/c2 0. Table HTl compares our measurement of 
Am with the available results from the MARK I, DELCO 
and MARK II Collaborations , p^S] ; our measurements 
agree with the MARK I and DELCO results and are 
~ 3(7 above the MARK II result. The systematic error 
on the mass measurement is evaluated by varying the 
background function, the width of the 0(3770) within its 
errors 0, changing the fit range and changing the bin 
width. It also includes the uncertainty in the mass. 

A fit to the M{D+D-) distribution of Fig.j^b) yields 
7.7 ± 4.2 events. Here the 0(3770) mass was fixed at 
3778.4 MeV/c^, the value found from the M{D"D") fit. 

We studied a MC sample of generic BB events that 
has the same size as the data sample. We find that 
the DD invariant mass exhibits a smooth behavior with- 
out peaks. We also analyzed off-resonance data taken 
60 MeV below the T(45) with a_ 10 fb-^ data sample. 
The same selection applied to D^D^K^ combinations re- 
sults in one event over the whole M{D^D^) region, which 
corresponds to a negligible contribution from the contin- 
uum. 

The 0(3770) D°D° helicity distribution, deter- 
mined by fitting the M{D^D^) distribution for the 
0(3770) yield in each of eight cos0^(377o) bins 0, is 
shown in Fig. 0| The points are data and the histogram 
gives the result of a fit using MC-based expectations for 
a J^'-^ = 1 0(3770) with fioated normalization. The 
confidence level of the fit is 10.5%. 

The MC-determined efficiencies for B+ 0(377O)ii:+ 
followed by 0(3770) D^D^ and V(3770) ^ D+ D- 
are 10.3 x lO^'' and 5.7 x 10^**, respectively. This 
gives B{B+ V(3770)i^+) x 6(0(3770) ^ D"D") = 
(0.34 ± 0.08 ± 0.05) X 10"^ and B{B+ -> 0(377O)if +) x 




M(DD) (GeV/c^) 



FIG. 2: a): The M{D°D°) distribution for the events from 
the _B-signal region. The dashed line represents the back- 
ground parameterization (see the text). The hatched his- 
togram is constructed from the AE sidebands, b): Fitted 
M{D+D~) distribution. 



6(0(3770) D+D-) = (0.14±0.08±0.02)xl0-^ where 
the first error is statistical and the second is systematic. 
In addition to the sources already mentioned, the sys- 
tematic error comes from varying the signal and back- 
ground shapes in M {DD) fitting (5%) and from vary- 
ing the 7r° reconstruction efficiency (6%). From these 
two measurements we obtain the ratio 6(0(3770) 
D°D°)/B{ij{3770) D+D-) = 2.43±1.50±0.43. Given 
the large errors, our measurement is consistent with the 
previous measurement of this rati o by the MARK III 
Collaboration of 1.36 ± 0.23 ± 0.14 jij. 

To extract B{B^ ■0(3770)/^"'") from the measure- 
ments of B{B+ -> 0(377O)X+) X 6(0(3770) D°D°) 
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FIG. 3: Helicity distribution for B+ ^(3770)7^+ decay 
followed by ^(3770) D° D° . The points with errors are 
obtained from fits to the M(D'^D^) data in each cosS^(377o) 
bin. The histogram shows the expected distribution for 
il){ZnO)K+ (see text). 



and B{B+ -> ^{mQ)K+) x ^(^'(3770) ^ D+D-), we 
assume that the D'^D'^ and D^D^ modes totally saturate 
the '0(3770) decay width. Summing both measurements 
gives B{B+ iP{3770)K+) = (0.48±0.11±0.07) x lO^^. 

Belle recently reported the observation of a nar- 
row charmonium-like state X(3872) that decays to 
Tr^Tr~J/tp HJ. This state, which is seen in the ex- 
clusive decay B — > KX{3872), is above DD thresh- 
old. Information about the X(3872) DD decay rate 
would be useful for determining its quantum numbers. 
We refitted the D°D° and D+D' invariant mass dis- 
tributions including possible contributions from B~^ — > 
X{3872)K+, X(3872) ^ DD decays_. The fits yield 2.1± 
1.8 and 0.4 ± 0.8 events for the D°D° and D+D' chan- 
nels, respectively. From this we determine 90% CL upper 
limits B{B+ X{3872)K+) x S(X(3872) -> D'^D^) < 
6 X 10-5 and B{B+ X(3872)X+) x B{X{3872) -> 
D^D~) < 4 X 10^5 . We have also searched for a pos- 
sible reflection from B+ X{3872)K+, X{3872) 
D^D^TT^ decays. This decay mode of the X(3872) is 
interesting because it is predicted to be large if the 
X(3872) is a DD* multiquark "molecular state" 
A MC study shows that these decays produce a nar- 
row, nearly gaussian reflection peak {a — 9 MeV) cen- 
tered at AE = --145 MeV. Using the D°D°K+ signal 
described above, we require M{D^D^) to be less than 
Af(X(3872)) - M(7rO) = 3737 MeV/c^ and fit the re- 



sulting AE distribution to a gaussian with mean and 
width fixed at the values expected for the reflection peak 
and a linear background contribution. The fit yields 
2.2±1.7 events. From this we determine a 90% CL upper 
limit B{B+ X{3872)K+) x B{X{3872) D^D^tt^) < 
6 X IQ-^. 

In summary, we have measured the branching fraction 
for B+ D°D°K+ decay to be B{B+ D°D°K+) = 
(1.17±0.21±0.15) X 10-3. A search for B+ D+ D- K+ 
decay results in an upper limit oiB{B^ D^ D^ K^) < 
0.90 X 10-3 (90% CL). Wc observe a peak_in the D"D° 
invariant mass spectrum from B^ D'^D'^K^ decays 
with a mass near 3770 MeV/c^ that we attribute to ex- 
clusive i?+ ip{3770)K'^ decay. This signal, which 
has a statistical significance of 5.9 a, is the first obser- 
vation of this decay mode. The mass of the -0(3770) 
is measured to be 3778.4 ± 3.0 ± 0.8 MeV/c^ The 
value of B{B+ V(3770)if+) x S(0(377O) ^ D^D'^) 
is measured to be (0.34 ± 0.08 ± 0.05) x IQ-^. For 
B+ 0(377O)A"+ followed by 0(3770) ^ D+D- we 
extract B{B+ 0(377O)if +) x S(0(377O) -> D+D-) = 
(0.14 ± 0.08 ± 0.02) X 10-3. The ratio sltS^fo-} 
is 2.43 ± 1.50 ± 0.43. By assuming that the D°D° and 
D^D^ modes totally saturate the 0(3770) decay width 
we obtain ;B(B+ ^ V(3770)i4:+) = (0.48 ± 0.11 ± 0.07}x 
10-3 which is comparable to B{B+ 0(25')if+) ^. 
This result suggests a large amount of S-D mixing in the 
0(3770). 

For the decays B+ X{3872)K+ followed by 

X(3872) -> D°D° and D+D- we have set 90% CL upper 
limits on B{B+ X{m2)K+) x S(X(3872) ^ DD) 
of 6 X 10-5 and 4 x 10"^ respectively. For the decay 
B+ X{3872)K+ followed by X(3872) ^ D^D^tt^ we 
have set a 90% CL upper hmit of 6 x 10-^. 
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